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Treating the cosmological constant as a thermodynamic pressure and its conjugate quantity as 
a thermodynamic volume, we reconsider the critical behaviour of charged AdS black holes. We 
complete the analogy of this system with the liquid-gas system and study its critical point, which 
occurs at the point of divergence of specific heat at constant pressure. We calculate the critical 
exponents and show that they coincide with those of the Van der Waals system. 
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I. INTRODUCTION 



Thermodynamic properties of black holes have been 
studied for many years. During that time it turned out 
that black hole spacetimes can not only be assigned stan- 
dard thermodynamic variables such as temperature or 
entropy, but were also shown to possess rich phase struc- 
tures and admit critical phenomena, in complete analogy 
with known non-gravitational thermodynamic systems 
elsewhere in nature. Of especial interest is black hole 
thermodynamics in the presence of a negative cosmolog- 
ical constant. Asymptotically AdS black hole spacetimes 
admit a gauge duality description and are described by 
dual thermal field theory. Correspondingly one has a 
microscopic description of the underlying degrees of free- 
dom at hand. This duality has been recently exploited 
to study the behaviour of quark-gluon plasmas and for 
the qualitative description of various condensed matter 
phenomena. 

The history of the study of thermodynamic proper- 
ties of AdS black holes began with the seminal paper of 
Hawking and Page [1| , who demonstrated the existence of 
a certain phase transition in the phase space of the (non- 
rotating uncharged) Schwarzschild-AdS black hole. Since 
then our understanding of phase transitions and criti- 
cal phenomena has been extended to a variety of more 
complicated backgrounds, e.g., @, Q. Of particular in- 
terest is the discovery of the first order phase transition 
in the charged (non-rotating) Reissner-Nordstrom-AdS 
(RN-AdS) black hole spacetime 0, Q. This transition 
displays classical critical behaviour and is superficially 
analogous to a liquid-gas phase transition ala Van der 
Waals. We critically review this analogy in the appendix. 

Recently, the idea of including the variation of the cos- 
mological constant A in the first law of black hole thermo- 
dynamics has attained increasing attention From 



a general relativistic point of view, such a variation is a 
slightly awkward thing to do, as the cosmological con- 
stant should be treated as a fixed external parameter of 
the theory. Moreover, varying A in the first law corre- 
sponds to comparing black hole ensembles with different 
asymptotics. This is different from standard thermody- 
namic considerations, where the black hole parameters 
are varied in a 'fixed AdS background'. However, there 
exist at least three good reasons why the variation of 
A should be included in thermodynamic considerations: 
i) One may consider 'more fundamental' theories, where 
physical constants, such as Yukawa couplings, gauge cou- 
pling constants, Newton's constant, or the cosmological 
constant are not fixed a priori, but arise as vacuum ex- 
pectation values and hence can vary. In such a case it is 
natural to include variations of these 'constants' in the 
first law, e.g., u ) More pragmatically, in the pres- 

ence of a cosmological constant the first law of black hole 
thermodynamics becomes inconsistent with the Smarr re- 
lation (the scaling argument is no longer valid) unless the 
variation of A is included in the first law Q • A similar sit- 
uation occurs, for example, for Born-Infeld black holes, 
where the variation of the maximal field strength b must 
be included in the first law for it to be consistent with 
the corresponding Smarr relation [TB - [l7| . iii) Once the 
variation of A is included in the first law, the black hole 
mass M is identified with enthalpy rather than internal 
energy [7] . It is also natural to consider a variable ther- 
modynamically conjugate to A. Since A corresponds to 
pressure, its conjugate variable has dimensions of vol- 
ume, characteristic for a given black hole spacetime. In 
particular, using geometric units Gjv = fr = c = k = 1, 
in the case of an asymptotically AdS black hole in four 
dimensions, one identifies the pressure with 



P 



(1.1) 
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1 In fact, this idea was first applied (and is now commonly used) in Q. 
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um€ with V = (j^p) g q j ■ In the simplest case of a 
Reissner-Nordstrom-AdS (RN-AdS) black hole the vol- 
ume is given by 



V 



+ > 



(1.2) 



with r + being the radius of the black hole event horizon 
expressed in terms of the 'Schwarzschild radial coordi- 
nate'. More generally, the thermodynamic volume for a 
wide variety of known black hole solutions was studied in 
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It turns out that this volume seems to possess some 
universal properties and is conjectured to satisfy the re- 
verse isoperimetric inequality. If true, such an inequality 
may restrict some processes in black hole spacetimes. 

Taking these considerations seriously, once the black 
hole pressure and volume are identified, one may proceed 
to calculate various thermodynamic quantities employ- 
ing standard thermodynamic machinery. For example, 
one may study adiabatic compressibility, specific heat at 
constant pressure, or even the 'speed of sound' associ- 
ated with the black hole [§4l0j. Moreover, as noted by 
Dolan Q , this also opens an interesting possibility of re- 
considering the critical behaviour of AdS black holes in 
an extended phase space, including pressure and volume 
as thermodynamic variables. In particular, Dolan stud- 
ied P = P(V, T) equation of state for a rotating charged 
AdS black hole, observed an analogy with the Van der 
Waals P — V diagram, and determined its critical point. 

In this paper we explore this issue further and elabo- 
rate on Dolan's results. For simplicity, we concentrate on 
studying P — V criticality of (spherical) RN-AdS black 
holes. Our aim is to complete the identification with the 
liquid-gas system in the extended phase space, and to 
clarify the relationship with the previously observed Van 
der Waals-like behaviour of the system discussed in 0, Q . 
For this purpose, we study the behaviour of the Gibbs 
free energy in the fixed charge ensemble. Equipped with 
this we obtain the phase transition in the (P, T)-plane, 
standard for the liquid-gas transition. While this be- 
haviour is very similar to that uncovered in 0, H[, the 
results we obtain have a rather different interpretation 
~ they are not an analogy, but actually compare the 
'same physical quantities'. For example we find, as in 
the liquid-gas system, that the transition occurs in the 
(P, T)-plane and not in the (/3, Q)-plane. We also study 
the critical exponents and show that they coincide with 
those of the Van der Waals system. This completes the 
identification of the charged black hole first order transi- 
tion with the standard liquid-gas phase transition. 

Our paper is organized as follows. In the next section 
we review some basic facts about the Van der Waals fluid 
and its critical exponents. In Sec. 3 we demonstrate 
analogous behaviour of the RN-AdS black hole system 
to the Van der Waals fluid. Specifically, P — V critical- 
ity, Gibbs free energy, first order phase transition, and 
the behaviour near the critical point are identified with 
the liquid-gas system. Sec. 4 is devoted to conclusions. 
To compare our results with the Van der Waals analogy 
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FIG. 1. P - V diagram of Van der Waals fluid. The 

temperature of isotherms decreases from top to bottom. The 
two upper dashed lines correspond to the "ideal gas" phase 
for T > T c , the critical isotherm T = T c is denoted by the 
thick solid line, lower solid lines correspond to temperatures 
smaller than the critical temperature; T — To isotherm is also 
displayed. The constants a and b in Eq. (|2.1[) were set equal 
to one. 



studied in 0, Q we include an appendix. 



II. VAN DER WAALS FLUID 

Van der Waals equation, e.g. [18| . is a popular closed 
form modification of the ideal gas law, which approxi- 
mates the behaviour of real fluids, taking into account 
the nonzero size of molecules and the attraction between 
them. It is often used to describe basic qualitative fea- 
tures of the liquid-gas phase transition. The equation 
reads 



(p +£)(.-») 



kT. 



(2.1) 



Here, v = V/N is the specific volume of the fluid, P its 
pressure, T its temperature, and k is the Boltzmann con- 
stant. The constant b > takes into account the nonzero 
size of the molecules of a given fluid, whereas the con- 
stant a > is a measure of the attraction between them. 
Equivalently, one can expand this equation to write it as 
a cubic equation for v, 

Pv 3 - {kT + bP)v 2 + . 



ab = 0. 



(2.2) 



The qualitative behaviour of isotherms in the P — V 
diagram is depicted in Fig. 1. The critical point occurs 
when P = P(v) has an inflection point, i.e., 
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FIG. 2. Maxwell's equal area law. The 'oscillating' 
(dashed) part of the isotherm T < T c is replaced by an isobar, 
such that the areas above and below the isobar are equal one 
another. 



at the critical isotherm T = T c Writing (f272l) in the 
form P c (v — v c ) 3 = , and comparing the coefficients we 
conclude that 
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We also note that 



kT c = , v c — 3b , P c 
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(2.4) 
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is a universal number predicted for all fluids (independent 
of constants a and b). Defining 



P 



T 



Pc Vc ±c 



(2.6) 



we obtain the so called law of corresponding states 



8r = (3u-l)[p+ — 



(2.7) 



This law is universal and is valid under more general 
assumptions than the original derivation of the Van der 
Waals equation. 

For T < T c there is a liquid-gas phase transition in 
the system. To describe this phase transition one has to 



2 Note that for T < T 



3¥§ = §5^°' a P ar t °f the isotherm 
corresponds to a negative pressure. Similar to the black hole 
case discussed in the next section, this is unphysical and once 
the Maxwell's equal area law, 112.81 1. is employed this is replaced 
by the corresponding isobar of positive pressure. 



FIG. 3. Gibbs free energy of Van der Waals fluid. This 
picture shows the characteristic swallowtail behaviour of the 
Gibbs free energy as a function of pressure and temperature. 
This corresponds to a first-order liquid-gas phase transition 
which occurs at the intersection of G surfaces. The corre- 
sponding curve is called the coexistence line. We have set 
$ = 1. 



replace the 'oscillating part' of the isotherm by an isobar, 
according to Maxwell's equal area law 



vdP = . 



(2.8) 



as shown in Fig. 2. This prescription reflects the fact 
that at the transition both phases have the same Gibbs 
free energy, see (|2.9I) . Eq. (I2.8[) is an effective tool to find 
the pressure at which the phase transition occurs, or to 
calculate the change of the volume and the latent heat. 

To get more information about the phase transition 
we study the (specific) Gibbs free energy, G = G(P,T). 
For a fixed number of particles, this can be obtained by 
integrating its differential 



dG 



-SdT + vdP, 



(2.9) 



while using the Van der Waals equation. The unspecified 
integration function can be determined from comparing 
with (statistical) G of the ideal gas. The recovered spe- 
cific Gibbs free energy reads 



G = G(T, P) = -kT 1 + In 



(v - b)T 3 / 2 



h Pv , 



(2-10) 

Here, v is understood as a function of P and T, through 
the Van der Waals equation (|2.1[) . and $ is a (dimen- 
sionful) constant characterizing the gas. The Gibbs free 
energy is depicted in Fig. 3 and Fig. 4. 

The coexistence line of two phases, along which these 
two phases are in equilibrium, occurs where two surfaces 
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FIG. 4. Gibbs free energy of Van der Waals fluid. This 
figure depicts the qualitative behaviour of the Gibbs free en- 
ergy as a function of temperature for various pressures. The 
pressure decreases from right to left. The dashed line corre- 
sponds to P = 2P C , the thick solid line to P — P c , and the 
remaining solid lines display P — 0.6P C and P — 0.09P C , re- 
spectively. For P < P c there is a first-order transition in the 
system. 



of G cross each other. This line is depicted in Fig. 5. It 
is governed by the Clausius— Clapeyron equation 



dP 



coexistence 



S g — S[ 
v„ - vi 



(2.11) 



where S g , Si, and v g , uj, stand for the specific entropy 
and specific volume of the gas, liquid, phase, respectively. 
Alternatively, this line can be constructed by exploiting 
Maxwell's law (|2.8|l . or by finding a curve in the (P, T)- 
plane for which the Gibbs free energy and the Van der 
Waals temperature coincide for two different volumes vi 
and v g . 

Critical exponents describe the behaviour of physical 
quantities near the critical point. It is believed that they 
are universal, i.e., they do not depend on the details of 
the physical system^ Let us calculate a few of these 
exponents as predicted by the equation of corresponding 
states. We denote 



T -T c 

t = = T - 1 . 



= v-l. (2.12) 



The critical exponents a, ft, 7, 8 are defined as follows: 



3 They may, however, depend on the dimension of the system or 
the range of the interactions. In D > 4 dimensions they can be 
calculated by using the mean field (Landau's) theory, e.g., |18|| . 



FIG. 5. Coexistence line of liquid—gas phase. Fig. dis- 
plays the coexistence line of liquid and gas phases of the Van 
der Waals fluid in (P, T)-plane. The critical point is high- 
lighted by a small circle at the end of the coexistence curve. 



• Exponent a governs the behaviour of the specific 
heat at constant volume, 



C — T— 



(2.13) 



• Exponent j3 describes the behaviour of the order 
parameter 7/ = v g — vi (the difference of the volume 
of the gas v g phase and the volume of the liquid 
phase vi) on the given isotherm 



n = v g -V[(x \tf . 



(2.14) 



• Exponent 7 determines the behaviour of the 
isothermal compressibility Kt, 



Kt 



1 dv 

vdP 



<x\t\- 



(2.15) 



• Exponent <5 governs the following behaviour on the 
critical isotherm T = T c : 



IP -P c | OC \V-V C 



(2.16) 



Let us now, following [18], calculate these exponents 
for our system. To calculate C v we consider the specific 
free energy 



F(T, v) = G-Pv = -kT 1 + In 



{v - 6)T 3 / 2 



From here we can calculate the specific entropy 



v 

(2.17) 



(2.18) 
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and hence 



Os 

Cv=T df 



3k 

Y 



(2.19) 



Since this expression is independent of t we conclude that 
exponent a — 0. 

Expanding the equation of corresponding states, (|2.7|) . 
in the vicinity of the critical point we have 



p = 1 + it - 6i0 - -cp 3 + 0(tcj) 2 ,<j) 4 ) . 



(2.20) 



[The fact that we can neglect the term 0(t<p 2 ) in this ex- 
pansion is justified by Eq. (|2.23|) below.] Differentiating 
this for fixed t < we get 



dP 



(2.21) 



Let us denote <p g = (v g — v c )/v c the 'volume' of gas and 
similarly <pi the volume of the liquid. Using Maxwell's 
area law (|2 .8[) , and the fact that during the phase transi- 
tion the pressure remains constant, we have the following 
two equations: 



3 

p = 1 + it - 6t<f>t <pf = 1 + it - 6t<j) g 



= 



(2.22) 



These equations have a unique non-trivial solution, given 
by (j) g = —(pi — 2-\/— t. Hence we have 

V = v c {<Pg - (pi) = 2v c (j)g = iv c yf~t ji = 1/2 . 

(2.23) 

To calculate exponent 7, wc differentiate Eq. (|2.20j) to 



6P C t yv ' 



get 

dv 
dP 

Hence we have 

1 dv 
v oP 

We also note that using the formula 
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7 =1. 



Cp — C v 



-T 



dP 
df 



dv 



oc 



(2.24) 



(2.25) 



(2.26) 



That is, C p diverges at the critical point with the same 
critical exponent 7 as kt- 

Finally, the 'shape' of the critical isotherm t — is 
given by (j2~2T)]) ,i.e., 



p-1 



6 = 3. 



(2.27) 



This completes the calculation of basic critical exponents 
for the Van der Waals fluid. 



III. PHASE TRANSITION IN CHARGED ADS 
BLACK HOLE SPACETIME 

A. Charged AdS black hole 

To start with we review some basic thermodynamic 
properties of the spherical RN-AdS black hole. In 
Schwarzschild-like coordinates the metric and the U(l) 
field read 



ds 2 = -Vdt 2 
F = dA 



dP_ 
A=-Qdt. 



l dni 



(3.1) 
(3.2) 



Here, d£l 2 stands for the standard element on S 2 and 
function V is given by 



V = 1 



2M 
r 



Q 2 



(3.3) 



With this choice we have a solution of the Einstein- 
Maxwell system of equations following from the bulk ac- 
tion 



EM 



1 
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g [R — F 



(3.4) 



The position of the black hole event horizon is deter- 
mined as a larger root of V{r + ) = 0. The parameter M 
represents the ADM mass of the black hole and in our 
set up it is associated with the enthalpy of the system. Q 
represents the total charge. Using the 'Euclidean trick', 
one can identify the black hole temperature 



l = V'(r + ) 



1 



4-7T 47TJV 

and the corresponding entropy 



3r 2 
/2 



1 + ^-3- • ( 3 - 5 ) 



A = iirr 



+ ■ 



(3.6) 



The electric potential measured at infinity with re- 
spect to the horizon, is 



(3.7) 



The thermodynamic volume V of the black hole is given 
by (|1.2p and the corresponding pressure P by (jl.ip . With 
these relations, the solution obeys the first law of black 
hole thermodynamics in an extended (including P and V 
variables) phase space 

dM = TdS + $rfQ + VdP . (3.8) 

The corresponding Smarr relation 

M = 2(TS - VP) + $Q (3.9) 

can be derived from it by a scaling (dimensional) argu- 
ment Q. When P is treated as constant (i.e. if the 
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cosmological constant is not allowed to vary), (|3.8[) re- 
duces to the standard first law in the 'non-extended' 
phase space. In this case, however, the Smarr relation 
remains unchanged and no longer follows from the first 
law by the scaling argument. 

Previous work on the critical behaviour of RN-AdS 
black hole in the non-extended phase space demonstrated 
that in the canonical (fixed charge) ensemble, for Q < 
Q c and f3 z < (3 < f3 c , there exists a first order phase 
transition in the system (HQ. This occurs in the (/3, Q)- 
plane and the phase transition is quite analogous to the 
liquid-gas phase transition occurring in the (P, T)-plane 
for T < T c . We review this analogy in the appendix and 
find that it is not fully consistent. 

In what follows we concentrate on analyzing the phase 
transition of the AdS charged black hole system in the 
extended phase space while we treat the black hole charge 
Q as a fixed external parameter, not a thermodynamic 
variable. We shall find that an even more remarkable 
coincidence with the Van der Waals fluid is realized in 
this case. 



B. Equation of state 

For a fixed charge Q, Eq. (|3.5[) translates into the equa- 
tion of state for a charged AdS black hole, P = P(V, T), 



P 



T 



1 Q 2 

+ ~ — - , r + 



2r+ 8nr 2 + &irr\ 



W 

4-7T 



1/3 



(3.10) 



Here, V is the thermodynamic volume, given in terms of 
the event horizon radius r+, T is the black hole temper- 
ature, and Q its charge. 

Before we proceed further we perform dimensional 
analysis, translating the "geometric" equation of state 
(|3.10[) to a physical one. The physical pressure and tem- 
perature are given by 



Press = ^-P , Temp = '^-T , 
l 2 , 



he, 
T 



(3.11) 



while the Planck length reads Ip = ^f-. Multiplying 
(|3~TU|) with we get 



He he T 
Press = -pr-P — — 
l 2 P l P 2r+ 



fcTemp 
2l 2 P r+ 



(3.12) 



Comparing with the Van der Waals equation, (I2.1[) . we 
conclude that we should identify the specific volume v 
with 



v = 2l 2 P r+ 



(3.13) 



In other words it is the horizon radius r + , rather than 
the thermodynamic volume V, that should be associated 
with the fluid volume. 
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FIG. 6. P — V diagram of charged AdS black holes. 

The temperature of isotherms decreases from top to bottom. 
The two upper dashed lines correspond to the "ideal gas" one- 
phase behaviour for T > T c , the critical isotherm T — T c is 
denoted by the thick solid line, lower solid lines correspond to 
temperatures smaller than the critical temperature, and the 
T = To isotherm is also displayed. We have set Q = 1. 



Pursuing this identification further and returning to 
geometric units, the equation of state (|3.10[) can now be 
written at@ 



P= T --^ 
v 2ttv 2 



2Q2 rim 

— T j (3-15) 



The corresponding a P—V diagram" is depicted in Fig 6. 
Evidently, for Q ^ and for T <T C there is an inflection 
point and the behaviour is reminiscent of the Van der 
Waals gas. The critical point is obtained from 



dv 



d 2 P 

dv 2 



0, 



which leads to 
V6 



18ttQ 



2V6Q , P c 



96ttQ 2 



(3.16) 



(3.17) 



The critical radius corresponds to the critical thermody- 
namic volume 



8V6irQ 3 , 



(3.18) 



4 It is obvious that the criticality cannot happen in the 'grand 
canonical' (fixed $) ensemble. Indeed, in that case the equation 
of state reduces to the quadratic equation 

Pv 2 - 8nTv + 4(l - $ 2 ) = 0, (3.14) 

which does not admit critical behaviour. For a detailed discus- 
sion of thermodynamics in fixed <t> ensemble see [l9ll . 
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which will be used later for calculating the critical expo- 
nents. 

Inspecting the critical values (|3.17|) we find an inter- 
esting relation 



(3.19) 



which is exactly the same as for the Van der Waals fluid 
and is a universal number predicted for any RN-AdS 
black hole with arbitrary charge. Moreover the critical 
values (|3.17l) can be written in the form (|2.4p . with 



3 , 2V6Q 

a = — , o = 

4tt ' 3 



(3.20) 



It is the presence of the charge Q that makes the effective 
volume smaller, b oc Q. Defining further 

P v T 



V 



(3.21) 



equation of state (|3.10p translates into 'the law of corre- 
sponding states' given by, cf. Eq. (|2.7p . 



8t = 3i/ [p 



(3.22) 



Similar to the Van der Waals equation (see Footnote 2) 
there exists a temperature To, given by 



To = 



V3 
18ttQ 



(3.23) 



below which the pressure P becomes negative for some 
r_|_. As with the fluid this behaviour is unphysical, as 
one has to replace the oscillatory part of the isotherm by 
an isobar, according to Maxwell's area law (|2.8[) . This 
is supported by studying the Gibbs free energy of the 
system in the next subsection. 

Finally, let us mention that the critical behaviour is 
present only for a black hole with horizon of spherical 
topology. For different topologies the equation of state 
becomes 



P 



T 

2^ 



Q 2 



8vrrj 



87T7-4 



(3.24) 



where k = is the toroidal (planar) case and k = — 1 
is the higher-genus (hyperbolic) case. Obviously, since 
all the terms on the r.h.s. are now non-negative, critical 
behaviour cannot occur. Consequently there is also no 
temperature Tq. 




FIG. 7. Gibbs free energy of charged AdS black holes. 

Figure shows characteristic swallowtail behaviour of the Gibbs 
free energy as a function of pressure and temperature for a 
fixed Q = 1. The first-order small-large black hole transition 
occurs at a curve which is an intersection of G surfaces. 
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C. Gibbs free energy 

In order to find the partition function of the system let 
us calculate its Euclidean action. For a fixed charge Q, 
one considers a surface integral 

I s = — - ( d 3 xVhK - — ( (Pxy/hn a F ab A b . 
8tt Jqm 47T J dM 

(3.25) 



FIG. 8. Gibbs free energy of charged AdS black hole. 

Fig. depicts G as a function of temperature for fixed Q = 
1, and increasing pressure P/P c = 0.1 (black line) top to 
P/P c = 1.6 (red line). The green line corresponds to critical 
pressure P = P c w 0.0033. Obviously, for T < T c w 0.043 
there is a (small black hole)- (large black hole) first-order phase 
transition. 
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The first term is the standard Gibbons-Hawking term 
while the latter term is needed to impose fixed Q as a 
boundary condition at infinity. The total action is then 
given by 



I 



EM 



L 



(3.26) 



where I em is given by (|3.4p , and I c represents the invari- 
ant counterterms needed to cure the infrared divergences 
I2l| . The total action was first calculated in SI 
and reads 



3l 2 Q 2 



(3.27) 



In the standard approach this action is associated with 
the free energy, stressing its dependence on extensive 
quantity Q. However, since we are considering an ex- 
tended phase space, and the action has been calculated 
for fixed A, following [9j we associate it with the Gibbs 
free energy for fixed charge, 



G = G{T,P) = \ 



8tt 3 3Q 2 
r + - -Pr\ + — 



(3.28) 



Here, r + is understood as a function of pressure and tem- 
perature, r + = r+(P,T), via equation of state (|3.10[) . 

The behaviour of G is depicted in Fig. 7 and Fig. 8. 
Since the G surface demonstrates the characteristic 'swal- 
low tail' behaviour, there is a first order transition in the 
system. The coexistence line in the (P, T)-plane can be 
obtained from Maxwell's equal area law, the Clausius- 
Clapeyron equation (|2.1ip . or simply by finding a curve 
in the (P, T)-plane for which the Gibbs free energy and 
temperature coincide for small (with r + = r s ) and large 
(r+ = r{) black holes. This is depicted in Fig. 9. 

We note that the coexistence line in the (P, T)-plane 
and the one of the Van der Waals fluid in the previous 
section look very much alike. The small-large black hole 
transition occurs for T < T c , for any non-trivial value of 
the charge Q. Furthermore we note some important dif- 
ferences between this diagram and an analogous diagram 
displayed in the appendix. There the transition occurs in 
(Q, /3)-plane, for Q < Q c and (3 < (3 C . Moreover there is 
a lower bound, Zorro's temperature Pz, below which the 
transition no longer occurs. This analogy also 'confuses' 
extensive and intensive thermodynamic quantities — see 
the appendix for more details. 



D. Behavior near critical point 

Let us now compute the critical exponents a, S 
for the black hole system, following the discussion in the 
previous section. We start with the behaviour of the 
specific heat at constant thermodynamic volume CV ■ For 
this we consider the free energy 



F(T, V) = G - PV = - 



r+ 



2nTr 



+ 



^ ) . (3.29) 



Critical Point 



0.003 



0.002 




0.02 0.03 0.04 0.05 



FIG. 9. Coexistence line of charged AdS black hole. 

Fig. displays the coexistence line of small-large black hole 
phase transition of the charged AdS black hole system in 
(P, T)-plane. The critical point is highlighted by a small circle 
at the end of the coexistence line. 



From here we calculate the entropy 



S(T,V) 



dF 
df 



■ + ' 



(3.30) 



which coincides with (I3.6[) . Since this is independent of 
T we find that Cy = and hence exponent a = 0. 
Expanding around a critical pointO 



t= 1 . 

Tr 



Vr. 



(3.31) 



we approximate the law of corresponding states (|3.22j) as 

4 



8 8 r! 

p = 1 + -t tU) - —oj 3 + 0{tuj 2 ,uj 4 ) 

y 3 9 81 v ' ' 



(3.32) 



[The truncation of the series is justified by formula (|3.35|) 
below.] Differentiating this series for a fixed t < we get 



4 

27" 



dP = — -P r (Qt + uj 2 ) du. 



(3.33) 



Hence, employing Maxwell's equal area law (|2.8p . while 
denoting u s and uji the 'volume' of small and large black 
holes, we get the following two equations: 



4 

81" 



o o o 4 q 

p = 1 + — t — —tuii — — ojf = 1 H — t tui s , 



0=1 u}(6t + u 2 )du. 



(3.34) 



5 Note that we could also write ui' = — 1 and the same results 

for the critical exponents would be obtained. 
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These have a unique non-trivial solution given by lo s 
—lui = 3\/—2t. Hence we have found that 

1 



77 = V c (ui - u.) = 2V c ui = 6V c V-2t /3 >; . 

(3.35) 

To calculate the exponent 7, we differentiate Eq. (I3.32|) 
to get 



&V 
dP 

Hence we have 

VdP 



Kt 



9Vcl 

SP C t 



9 11 

8PW 



7 



(3.36) 



(3.37) 



The critical point corresponds to the divergence of the 
specific heat at constant pressure 



C p C v 



-T 



dP 
df 



dV 



12ttQ 2 
t 



(3.38) 



That Cp but not Cy diverges at our critical point can 
be also directly seen from the definition 



r -r dS 
Cp - T df 

and the expression [9] 



Cv ~ T dT 



1 



TrQ 2 
S 



8PS 



We obtain Cy = 0, and 



Cp = 2S 



8PS 2 + S-irQ 2 
8PS 2 -S + 3irQ 2 



(3.39) 



(3.40) 



(3.41) 



where stability requires Cp > 0. The specific heat Cp 
becomes singular at 8PS 2 — S + 3itQ 2 = 0, and since 



S = 7rr+ this occurs when 



8iVrt 



3ttQ 2 = , 



i.e., exactly at the critical point. 

Finally, the 'shape of the critical isotherm' t 
given by ((332]) 



p-l 



81 



-oj 



(3.42) 
= is 

(3.43) 



This completes the calculation of basic critical exponents. 
Comparing with the previous section, where these were 
calculated for the Van der Waals fluid, we find that all of 
them coincide. These coefficients can also be compared 
with critical exponents in non-extended phase space, e.g., 
[22^ . Since critical exponents may depend on the dimen- 
sionality of the system, it would be interesting to repeat 
our calculation for charged BTZ black hole or higher- 
dimensional RN-AdS black holes. 



IV. CONCLUSIONS 

In this paper we have studied the thermodynamic be- 
haviour of a charged AdS black hole in an extended 
phase space — treating the cosmological constant and its 
conjugate quantity as thermodynamic variables associ- 
ated with the pressure and volume, respectively. For a 
fixed black hole charge, this identification allowed us to 
write the equation of state as P = P(V,T) and study 
its behaviour using the standard thermodynamic tech- 
niques. We have shown that the system admits a first- 
order small-large black hole phase transition which in 
many aspects resembles the liquid-gas change of phase 
occurring in fluids. Namely, for any non-trivial value of 
the charge, there exists a critical temperature T c below 
which this phase transition occurs. Similar to fluids, one 
can associate with it a coexistence line in (P, T)-plane 
which terminates at the critical point behind which it is 
no longer possible to distinguish the two phases. We have 
further studied the behaviour of certain physical quanti- 
ties near the critical point and calculated the correspond- 
ing critical exponents. These were shown to coincide with 
those of the Van der Waals fluid. 

In many aspects our considerations are similar to those 
performed in [U, [Bj] ; for convenience we reviewed them in 
the appendix. Therein the critical behaviour of charged 
AdS black holes was studied in the non-extended phase 
space and was found similar to that of the Van de Waals 
fluid. However the similarity found in these papers is a 
mathematical analogy rather than an exact correspon- 
dence. We point out in the appendix that one has to, for 
example, identify the fluid temperature with the charge 
of the black hole, the volume with the horizon radius and 
so on, confusing intensive and extensive quantities. On 
the other hand, in our approach to the problem we actu- 
ally compare the same physical entities: the temperature 
of the fluid is identified with the Hawking temperature 
of the black hole, the volume of the fluid associates with 
the thermodynamic volume of the black hole etc. Conse- 
quently the coincidence with the fluid behaviour is more 
direct and precise. We find a very similar behaviour of 
isotherms in P — V plane, coexistence line in P — T plane, 
as well the same critical exponents. All these results seem 
to justify and support the idea that considering the ex- 
tended phase space, and hence treating the cosmological 
constant as a dynamical quantity, is a very interesting 
theoretical possibility. 

We close by noting that we have limited ourselves 
to studying only the very basic properties of the sim- 
plest possible four-dimensional non-rotating charged AdS 
black hole. We have not discussed stability issues, the 
calculation of fluctuations, or the construction of Lan- 
dau's theory as done in [5j. Also, we have left out con- 
sidering the impact of possible non-linear electrodynam- 
ics extensions of RN-AdS black hole, e.g., (23l - [25| . or the 
impact of rotation which, especially in higher dimensions 
where one has more rotation parameters, may make the 
presence of possible phase transition(s) more interesting. 
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Also interesting would be to pursue the study 26] , or to 
investigate the recently demonstrated isomorphism be- 
tween phase structure of hairy black holes in massive 
gravity and RN-AdS black hole (27[. All these are left 
as future possible directions. 
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Appendix: Van der Waals analogy 

The thermodynamics of charged AdS black holes in 
the non-extended phase space was first studied in [3, Q . 
The authors considered two thermodynamic ensembles, 
one with a fixed electric potential $ at infinity, the other 
with the fixed black hole charge Q at infinity. To capture 
the fact that $ is an intensive quantity, whereas Q is ex- 
tensive, they called the corresponding ensembles as grand 
canonical and canonical, respectively. It turned out that 
the canonical ensemble admits a first order small-large 
black hole phase transition which is in many respects 
analogous to the liquid-gas phase transition in fluids. 
This triggered wide interest and various comparisons to 
the Van der Waals critical behaviour were studied, e.g., 
[22l [281-30!] and references therein. In this appendix we 
recapitulate this interesting analogy and relate it to the 
results in the extended phase space performed in this pa- 
per. 

We shall regard the cosmological constant A, and hence 
thermodynamic pressure, as a constant and work in the 
"canonical ensemble", i.e., the action is calculated for 
a fixed charge Q. The analysis again starts with the 
equation of state (|3.5p , which is now written as 



P = p(r+,Q) = 



47T7\i 



1 



(A.l) 



or, using the fact that the 'electric potential' $ reads 
* = Q/r+, 



P = P($,Q) 



Z2$2_/2$4 +3 Q2 



(A.2) 



The Euclidean action of the canonical ensemble is as- 
sociated with the quantity G, (|3.28[) . To stress its depen- 
dence on the extensive quantity Q, we shall denote it in 
this appendix as F — F(Q,f3). This thermodynamic po- 
tential governs the behaviour of the canonical ensemble. 
It reads 



1 



F = F(Q,P) = — [l 



3l 2 Q 2 



(A.3) 



1 fl 2 Q Q 



4.1- 



7H (^-f3+ 3 W) - ( A - 4 ) 




F 0.2- 



FIG. 10. Thermodynamic potential F. The thermody- 
namic potential F is depicted as a function of charge and 
inverse temperature for fixed / = V3. The characteristic 
swallowtail behaviour corresponds to the presence of the first- 
order phase transition in the canonical ensemble. Such fea- 
tures are absent in the grand canonical (fixed <£>) ensemble. 



where r + and $ are understood as functions of Q and 
P, as given by Eqs. (jA.lj) and (|A.2|) . Surfaces of F 
are displayed in Fig. 10. Obviously, F demonstrates the 
characteristic swallowtail behaviour, which corresponds 
to the first-order phase transition in the canonical en- 
semble. [Similar behaviour, on the other hand, is absent 
in the grand canonical (fixed $) ensemble.] 

To obtain an analogy with the liquid-gas phase tran- 
sition in fluids, we have to identify F with the Gibbs free 
energy G — G(P, T) of the fluid. For this we have two 
options, both of which were suggested in [J, |5j. One pos- 
sibility is to identify (3 with the pressure of the fluid and 
Q with its temperature. We call this analogy 1: 



Analogy 1 


fluid 


AdS black hole 


temperature 
pressure 
volume 


Q 

P 

r+ 



(A.5) 



The other possibility is to identify P with the fluid tem- 
perature and Q with its pressure — analogy 2: 



Analogy 2 


fluid 


AdS black hole 


temperature 
pressure 
volume 


P 

Q 



(A.6) 



Let us study both these possibilities in greater detail. 
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Q 0.6 



FIG. 11. P — r+ diagram. This is an analogue of the p — V 
diagram of the fluid. Isocharge lines are displayed for various 
charges. The upper dashed line corresponds to the "ideal gas 
phase" for Q > Q c , the critical isocharge line Q — Q c is high- 
lighted by the thick solid line, lower solid lines correspond to 
Q < Qc for which the phase transition occurs. Cosmological 
constant was set I = 



FIG. 12. "Zorro" diagram. The thermodynamic potential 
F is depicted as a function of Q; I = y/3. Value of the inverse 
temperature P decreases from top to bottom. Upper dashed 
line corresponds to P > p c , the upper thick solid line to/3 = 
Pc the lower thick solid line to "Zorro's signature" [5[ at 
P — Pz- The phase transition occurs for Pz < P < Pc- The 
bottom dashed line displays the disconnected line of P < Pz- 



1. Analogy 1 

In this analogy, one identifies the fluid temperature 
with the black hole charge Q, and the fluid pressure with 
the black hole inverse temperature /?. Since the conjugate 
quantity to P is the black hole entropy S = nr\, which 
depends entirely on the horizon radius r + , it is natural to 
identify r + with the volume of the fluid, see table (IA.5I) . 
The equation of state (|A.ip . can now be written as the 
following quartic for v = r+: 



Pv 3 



4tt/ 2 



pi 2 



3w 



= 0. 



(A.7) 



which is going to be compared with (|2.2p . The corre- 
sponding P — r + diagram is depicted on Fig. 11. This is 
an analogue of p — V diagram of the fluid displayed in 
Fig. 1. 

The phase transition occurs for Q < Q c and Pz < P < 
P c . The coexistence line of the equilibrium between small 
and large black hole phase is depicted on Fig. 13. The 
critical point is given by 



dp 
dr+ 



= 0. 



d 2 p 

dr 2 . 



= 0, 



which happens for 

Qc = L 



I 



(A. 



(A.9) 



Contrary to the fluid (or extended phase space) case, 
there exists a temperature, "Zorro's" temperature pz, 



below which the phase transition no longer occurs. This 
temperature is the temperature of the last connected P = 
const curve displayed in Fig. 12. It is given by 



PZ = TTl. 



(A.10) 



Another point worth noting is the fact that the diagram 
F = F(Q), displayed in Fig. 12, significantly differs from 
the corresponding fluid diagram in Fig. 4. We see that 
the analogy is broken here. The reason is that the swal- 
lowtail of the potential F, displayed on Fig. 10, 'faces' the 
/3-axis rather than the Q-axis. Consequently in Fig. 13, 
where the coexistence line is displayed, the axes p and 
Q 'swap' their roles to what one would expect from the 
corresponding fluid Fig. 5. 



2. Analogy 2 

In this second analogy, one identifies the fluid pressure 
with the black hole charge Q, and the fluid temperature 
with the black hole inverse temperature p. Since the 
conjugate quantity to Q is the electric potential $, it is 
reasonable to identify it with the volume of the fluid, see 
table (|A.6|> . The equation of state (|A.2|) . can now be 
written as the following quartic for v = $: 

3Q3 
vl 2 



Qv 3 - Qi 



4ttQ 2 
PI 2 



= 0. 



(All) 



which is to be compared with (|2.2[) . The corresponding 
Q — $ diagram is depicted on Fig. 14. This is an analogue 
of the P — V diagram in Fig. 1. 
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0.4- Q 



0.3- 



Critical Point 



0.4- 



P>Pc 



0.2- 



0.1- 



LARGEBH SMALL BH 



0.2 0.4 0.6 



FIG. 13. Coexistence line in the (f3, Q)-plane. Here we 
illustrate the small-large black hole coexistence line in (/3, Q) 
plane for fixed I = \/3, analogous to the liquid-gas coexistence 
line in (p,T)-plane. 



The critical point occurs at 

I 1 ttIVG , k s 

Q c = - 6 , * c = ^, Pc = —- (A.12) 

The corresponding coexistence line is depicted again on 
Fig. 13. The difference is that now the axes have the 
"correct position". Also Fig. 15 is now "quite analogous" 
to Fig. 4. 



3. Summary 

Both outlined analogies have many interesting features 
which resemble aspects of liquid-gas phase transition in 
the fluid system. However, these analogies are not "ex- 
act" : there are certain features, such as the existence of 
Zorro's temperature, which are not present in the fluid 
system. Moreover, the identification itself is a mathe- 
matical analogy, rather than an identification of simi- 
lar physical quantities. In fact, one even has to identify 
quantities of 'opposite' thermodynamic character: in the 
first analogy one identifies the extensive quantity Q with 
the temperature of the fluid, which is an intensive quan- 
tity. In the second analogy, one even associates intensive 
pressure with Q and extensive volume with $. These un- 
pleasant features disappear when one uses the extended 
phase space analysis performed in the main text. 

Let us finally mention that similar to the main text, 
one can study the behaviour of physical quantities near 
the critical point and calculate the corresponding critical 
exponents. This was done recently in [22j . where it was 
found that the analogous critical exponents coincide with 



FIG. 14. Q — $ diagram. This is a second analogue of 
p — V diagram of the fluid. The inverse temperature f3 of 
isoterms decreases from top to bottom. The upper dashed 
line corresponds to the "ideal gas phase" for j3 > /3 C , the 
critical isoterm ft = j3 c is highlighted by the thick solid line, 
two lower solid lines correspond to f3z < P < /3 C for which the 
phase transition occurs, the bottom dashed line displays an 
isoterm at /3 < /3z- Cosmological constant was set I = v3. 



F 

0.4- Q>Qc 




FIG. 15. Potential F as function of /3. The thermody- 
namic potential F is depicted as a function of temperature for 
snapshots of fixed charge and I — \/3. The value of charge de- 
creases from top to bottom. Upper dashed line corresponds 
to Q > Q c , thick solid line to Q = Q c , solid line with 'lit- 
tle bill' to < Q < Q c - The last solid line corresponds to 
the Schwarzschild-AdS case, Q = 0. The radius r+ of black 
holes increases from right to left. In the Schwarzchild case 
the 'branch of small black holes' disapears completely (it co- 
incides with ,0-axis and represents r+ = black holes). 
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those of the Van der Waals fluid and are independent of the number of dimensions of the RN-AdS black hole. 



S. Hawking and D. N. Page, Thermodynamics of Black [16 
Holes in anti-De Sitter Space, Commun.Math.Phys. 87 
(1983) 577. 

M. Cvetic and S. S. Gubser, Phases of R charged black [17 

holes, spinning branes and strongly coupled gauge 

theories, JHEP 9904 (1999) 024, [hep-th/9902195] . 

M. Cvetic and S. Gubser, Thermodynamic stability and [18 

phases of general spinning branes, JHEP 9907 (1999) 

010, hep-th/9903132] . 

A. Chamblin, R. Emparan, C. Johnson, and R. Myers, [19 
Charged AdS black holes and catastrophic holography, 
Phys.Rev. D60 (1999) 064018, [hep-th/9902l70] . 

A. Chamblin, R. Emparan, C. Johnson, and R. Myers, 
Holography, thermodynamics and fluctuations of charged [20 
AdS black holes, Phys.Rev. D60 (1999) 104026, 
[hep-th/9904197] . 

M. M. Caldarelli, G. Cognola, and D. Klemm, 
Thermodynamics of Kerr-Newman-AdS black holes and [21 
conformal field theories, Class. Quant. Grav. 17 (2000) 
399-420, hep-th/9908"022"] . [22 
D. Kastor, S. Ray, and J. Traschen, Enthalpy and the 
Mechanics of AdS Black Holes, Class. Quant. Grav. 26 
(2009) 195011, _ arXiv:0904.2765l . [23 

B. Dolan, The cosmological constant and the black hole 
equation of state, Class. Quant. Grav. 28 (2011) 125020, 
|arXiv:1008.5023] . 

B. P. Dolan, Pressure and volume in the first law of [24 
black hole thermodynamics, Class. Quant. Grav. 28 
(2011) 235017, [arXiv:1106.6260] . 

B. P. Dolan, Compressibility of rotating black holes, [25 

Phys.Rev. D84 (2011) 127503, [arXiv: 1109.0198] . 

M. Cvetic, G. Gibbons, D. Kubiznak, and C. Pope, 

Black Hole Enthalpy and an Entropy Inequality for the [26 

Thermodynamic Volume, Phys.Rev. D84 (2011) 024037, 

|arXiv: 1012.28881 . 

H. Lu, Y. Pang, C. N. Pope, and J. F. Vazquez-Poritz, [27 
AdS and Lifshitz Black Holes in Conformal and 
Einstein- Weyl Gravities, larXiv: 12 04. 1062 

G. W. Gibbons, R. Kallosh, and B. Kol, Moduli, scalar [28 

charges, and the first law of black hole thermodynamics, 

Phys.Rev. Lett. 77 (1996) 4992-4995, [hep-th/9607T08] . 

J. Creighton and R. B. Mann, Quasilocal [29 

thermodynamics of dilaton gravity coupled to gauge 

fields, Phys.Rev . D52 (1995) 4569-4587, 

[gr-qc/9505007] . 

D. Rasheed, Nonlinear electrodynamics: Zeroth and first [30 
laws of black hole mechanics, jhep-th/970208T| 



N. Breton, Smarr's formula for black holes with 
non-linear electrodynamics, Gen. Rel. Grav. 37 (2005) 
643-650, [gr-qc/0405116] . 

W. Yi-Huan, Energy and first law of thermodynamics 
for Bornlnfeld AdS black hole, Chinese Phys. B19 
(2010) 090404. 

N. Goldenfeld, Lectures on Phase Transitions and the 
Renormalization Group, (Westview Press, New York, 
1992). 

C. Peca and J. Lemos, Thermodynamics of 
Reissner- Nordstrom anti-de Sitter black holes in the 
grand canonical ensemble, Phys.Rev. D59 (1999) 
124007, [gr-qc/9805004] . 

R. Emparan, C. V. Johnson, and R. Myers, Surface 
terms as counterterms in the AdS / CFT 
correspondence, Phys.Rev. D60 (1999) 104001, 
hep-th/9903238 . 

R. B. Mann, Misner string entropy, Phys.Rev. D60 
(1999) 104047, [hep-th/9903229] . 

C. Niu, Y. Tian, and X. Wu, Critical Phenomena and 
Thermodynamic Geometry of RN-AdS Black Holes, 
larXiv: 1104.30661 

S. Fernando and D. Krug, Charged black hole solutions 
in Einstein- Born-Inf eld gravity with a cosmological 
constant, Gen. Re l. Grav. 35 (2003) 129-137, 
[hep-th/030"6120] . 

R. Banerjee and D. Roychowdhury, Critical phenomena 
in Born-Inf eld AdS black holes, Phys.Rev. D85 (2012) 
044040, [arXiv:1111.0i"47] . 

R. Banerjee and D. Roychowdhury, Critical behavior of 
Born Infeld AdS black holes in higher dimensions, 
larXiv:1203.0TT8l 

P. Nicolini and G. Torrieri, The Hawking-Page 
crossover in noncommutative anti-deSitter space, JHEP 
08 (2011) 097, llarXiv:1105.0188| . 
F. Capela and G. Nardini, Hairy Black Holes in 
Massive Gravity: Thermodynamics and Phase 
Structure, larXiv:1203.4222l 

A. Sahay, T. Sarkar, and G. Sengupta, Thermodynamic 
Geometry and Phase Transitions in Kerr-Newman-AdS 
Black Holes, JHEP 1004 (2010) 118, |arXiv:1002.2538] . 
A. Sahay, T. Sarkar, and G. Sengupta, On the 
Thermodynamic Geometry and Critical Phenomena of 
AdS Blac k Holes JHEP 1007 (2010) 082, 
[arXiv:1004.1625| . 

A. Sahay, T. Sarkar, and G. Sengupta, On The Phase 
Structure and Thermodynamic Geometry of R- Charged 
Black Holes, JHEP 1011 (2010) 125, |arXiv : 1009.2236 . 



